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ABSTRACT 

We investigate the Statistical Isotropy and Gaussianity of the CMB fluctuations, using 
a set of multipole vector functions capable of separating these two issues. In general 
a multipole is broken into a frame and 21 — 3 ordered invariants. The multipole frame 
is found to be suitably sensitive to galactic cuts. We then apply our method to real 
WMAP datasets; a coadded masked map, the Internal Linear Combinations map, and 
Wiener filtered and cleaned maps. Taken as a whole, multipoles in the range £ = 2 — 10 
or £ = 2 — 20 show consistency with statistical isotropy, as proved by the Kolmogorov 
test applied to the frame's Euler angles. This result in not inconsistent with previous 
claims for a preferred direction in the sky for £ = 2, ...5. The multipole invariants 
also show overall consistency with Gaussianity apart from a few anomalies of limited 
significance (98%), listed at the end of this paper. 

Key words: Cosmic microwave background - Gaussianity - Statistical Isotropy 



1 INTRODUCTION 

The recent WMAP data has enabled analysis of the CMB to 
much higher precision than ever before. The CMB acts like a 
curtain on the early universe and a back light on later times. 
It is also the largest thing we can observe, and therefore is a 
strong probe of features on cosmological scales. Inflationary 
theories predict the CMB fluctuations to be Gaussian 
and Statistically Isotropic (SI); the WMAP data should 
be the best way to test these predictions. Unfortunately, 
there ha ve been reports of quite severe anomalies in this 
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onset of a preferred direction. Such statistical anisotropy 
would be in conflict with inflationary theories, in the 
sense that any inflation model which left the universe 
just large enough to observe any anisotropy induced by, 
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for example, a non-trivial topology IWeeks fc al 120041: 
[ Tegmark. Oliveira-Costa. Hamiltonl 12003 : iRoukema&^Ll 
120041: ICornishfcal 12004jK would require extreme fine 
tuning. This aside, it is interesting to consider possible 
origins of anisotropy in the universe: p rimordial magnetic 
fields, strings or walls jBerera fc all20Q3lj . or an intrinsically 
inhomogeneous Universe jMoffatll2005ir 

In this paper we build on the wor k 

of iTegmark. Oliveira-Costa. Hamiltonl J2003J); 

ICopi. Huterer. Starkmanl J2003D : lLand. Magueiiol (l2004bl) . 
re-examining the issue of multipole alignment in the light of 
a formalism based on the Maxwell multipole vectors. These 
are an elegant way of extracting the extra degrees of freedom 
beside the pow er spectrum contained in a given multipole. 
As explained in lLand. Magueiiol l2004bl) . it is important to 
convert these vectors into 21? — 3 invariants (a generalisation 
of the concept of eigenvalue) and an orthonormal triad of 
vectors (a generalisation of the concept of eigenvectors of 
a hermitian matrix). Such an operation neatly separates 
the issues of non-Gaussianity and of violati ons of statistical 
isotropy . We proposed such a scheme in lLand. Magueijcl 
(2004b), and a further study appeared in iDennia (|2p04b), 
with analytical distributions provided. 

The plan of our paper is as follows. In Section we 
distinguish SI and non-Gaussianity in real and harmonic 
(Fourier) space. Then in Section[3]we review previous work 
done with multipole vectors, their invariants and frames, and 
present our new proposal. In Section 0] we demonstrate the 
effectiveness of the method with simulations, and we apply 
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these techniques to a variety of renditions of the WMAP 
data in Section |S| In the face of some concrete problems 
raised by real life data we conclude with an appraisal of 
our method, and how it stands in light of recent claims for 
anomalies in Section [(J T his paper is to be seen as a tech- 
nical companion paper to lLand. Magueiiol |^005) where the 
focus is more heavily placed on the analysis of the reported 
anomalies. 



2 STATISTICAL ISOTROPY AND 
GAUSSIANITY 

First we outline the definitions and implications of Statisti- 
cal Isotropy and Gaussianity in terms of the terminology of 
Real and Fourier space. 



2.1 Real space 

We consider the the temperature of the CMB photons re- 
ceived from different directions, ri, in terms of the fluctu- 
ation AT(n) = T(n) — T. Due to the random nature of 
the quantum fluctuations that seeded the initial perturba- 
tions (according to current standard cosmological models), 
our theories of this process can only determine the statis- 
tical behaviour of AT(n). A theory can only make definite 
predictions for ensemble averages: averages taken over an 
infinity of universes. 

One way to characterise all the statistical information 
of the AT field is through the n-point correlation functions: 



(AT(ni)AT(n 2 )...AT(n„)) 



(1) 



Note that without SI, these functions depend on 
rii,ri2, Tin- 

A function, f(r), is statistically isotropic if its statistical 
properties are invariant under translations r — > r + 5r , and 
there is no preferred direction trlaiian. SouradeerJ2005() . On 
the surface of a sphere, this is equivalent to the statistical 
properties being rotationally invariant. Statistical Isotropy 
of the temperature fluctuations therefore demands that the 
n-point correlation functions are invariant under rotations. 
They now only depend on the magnitude of the separation 
between the vectors, [r€i — n 3 * |. For example, the 2-point func- 
tion only depends on the distance between the 2 vectors, or 
equivalently the dot product: 



(AT(ni)AT(na)) = C(ni.n 2 ) 



(2) 



As mentioned, our models of the early universe can 
make predictions for the statistical properties of the field 
AT(h). The simplest models of inflation predict the AT(ii) 
to be a Gaussian random field. For a Gaussian field the n- 
point correlation functions are zero for odd n, and can be 
expressed in terms of the 2-point correlation function for 
even n. Therefore, the two-point correlation function con- 
tains all the statistical information for a Gaussian field. 

The combination of SI and Gaussianity means the two 
point correlation function contains all the information AND 
it is rotationally invariant. 



2.2 Fourier Space 

The usual approach to analysing CMB data, is to expand 
the function AT(h) in terms of the spherical harmonics, and 
therefore work with the spherical harmonic coefficients; ae m . 



AT(n) = ^2 a ^ 



Yim (n) 



(3) 



We can translate properties of the function AT(n) into 
properties of the ae m - As mentioned in Section 12.11 SI of 
AT(h) requires the 2-point correlation function of AT(n) 
to be rotationally invariant. The consequence for the aem 
is tha t the 2-point correlation function be of the form 
I2OOID : 

(ae imi ai 2nl2 ) = Ci 1 5t 1 i 2 8m 1 m 2 (4) 

We call Ce the power spectrum. 

Similarly, requiring the 3-point correlation function in 
real space to be SI d emands th at the 3-point function of the 
ai m s be of the form jHull200ll) : 



(^imi ^2 m 2 ^^3 m 3 ) 



mi 



V 2 

rn-z 



£3 
ni3 



B. 



(5) 



We call Bi 1 ( 2 £ 3 the Bi-spectrum. 

Gaussianity of the temperature fluctuations AT(h) de- 
mands that the coefficients, aim, are gaussian random vari- 
ables. Therefore, all the information about the field is con- 
tained in the 2-point correlation function {at imi at 2m3 }. All 
higher n-point correlations functions are to be obtained from 
the 2-point function by means of Wick's theorem. 

Therefore, in Fourier space terminology, the combina- 
tion of SI and Gaussianity means the power spectrum Ce 
contains all the information of the field. 



3 MULTIPOLES AND THE MULTIPOLE 
VECTORS 

A multipole, Te, represents all the fluctuations in AT of a 
certain scale, £ ~ 1/8. That is: 



AT(n) = 5^T 4 (ft) 

e 

Ti{n) = atmYtm{r 



(6) 
(7) 



The reality of AT means ae~ m is determined by ae m and 
Im(afo) = 0. There are therefore 2£ + 1 degrees of freedom 
in multipole Te. 

As discussed, SI demands that the n-point correlation 
functions be rotationally invariant, and in Fourier space 
this has consequences for the form the n-point functions 
(a^mi —■o-tn'mn) can take- In Equation^Jwe can see that SI 
demands that the at, m s depend only on I (no m-preference) 
and are independent for different £ (no inter-^ correlations). 
If we think about it in terms of the p.d.f of each ae m , the 
assumption of Gaussianity means that the shape of these 
p.d.f are Gaussian, and the assumption of SI means that 
their scaling are the same for the ae-mS of the same multi- 
pole, and independent for different multipoles. In fact, one 
can see that m-preference would indicate a preferred frame 
because the ae m s are not frame independent, they depend 
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on the frame of the Fourier Transform. Therefore, any m- 
preference would be related to a particular frame. 



3.1 Separation of degrees of freedom for the 
quadrupole 

A multipole is an irreducible representation of the SO (3) 
group. Therefore, from its 2£ + 1 degrees of freedom we 
should be able to extract 3 rotational degrees of freedom, 
i.e., an orthonormal frame. SI requires no inter-f correla- 
tions, and therefore the orthonormal frames of each multi- 
pole are independent. Equivalently, SI requires these frames 
to be uniformly distributed. The remaining 2£ — 2 degrees of 
freedom (d.o.f) will be invariants. Of these, one can be the 
overall power in the multipole Ci (a Gaussian d.o.f). The 
remaining 2^ — 3 invariants would probe m-preference and 
the Gaussianity of that multipole (these invariants should 
be independent for different multipoles). 

Therefore, if one separates the 3 rotational d.o.f from 
the invariants one can test SI in the form of no inter-£ 
correlations independently from the remaining issues of m- 
preference and Gaussianity. 

For the example of the quadrupole there is a simple way 
of separating the d.o.f jMagueiiolll99il : 

5T 2 = Qy*V (8) 

where Qij is a symmetric traceless matrix. From this ma- 
trix one may extract three eigenvectors and two indepen- 
dent combinations of invariant eigenvalues Ai (they sum to 
0). The eigenvalues d.o.fs are essentially the power spectrum 
C 2 (related to the sum of the squares of Ai) and the bispec- 
trum £2,2,2 (related to the determinant of the matrix, or 
the product of Ai). The eigenvectors provide us with the 
orthonormal frame. 



3.2 A general multipole 

In lLand. Magueiiol J2004bl) a new method was proposed 
that generalises this construction for any multipo le. It 
built on the work of lCopi. Huterer. StarkmarJ (120031) who 
first introduced the Multipole Vector formalism in CMB 
analysis. They showed that a multipole can be repre- 
sented in terms of I unit vectors and an overall magni- 
tude. The method has received much attention, from math- 
ema ticians securin g and exploring features of the formal- 
ism toenniJ2004alOlKatz. Week J200i lLaihieze-Re\l2004 
IWeeksll2004F ~ and from astrophysicists applying it to the 
WMA P data iCopi. Huterer. Starkmanll2003t ISchwarz fe"al 
2003). In fact t he formalism was first discovered by Maxwell 
( Maxwell 1891). He showed that for a real function f(x, y, z) 
which is an eigenfunction of the Laplacian on the unit sphere 
with eigenvalue — t{l + 1) (i.e., a spherical harmonic Y/ m ), 
there exists t unit vectors vi, Vf such that: 

f(x,y,z) = V Vl ...V Ve - (9) 

where V Vl = vi.V is the directional derivative operator, 
and r = \/ x 2 + y 2 + z 2 . A more us eful form of the repre- 
sentation is given by dDennisll2004aln 

/(r) = A(r.vi). .. (r.v £ ) +r 2 i? (10) 



where is A is a constant, and R is a complicated term fully 
defined by the components of the Vi. 

With the Multipole Vector method, the 2£ + l d.o.f are 
separated into I unit v ectors, Vi, and one invariant, A. In 
Land, Magueiiol J2004bh we proposed a further step of se- 
lecting 2 of the vectors as "anchor vectors" (we discuss this 
step further in Section f3. 311 . You then use these anchor vec- 
tors (L\,L 2 ) to define an orthonormal frame: 



ar=(Li-La) (11) 

y = (LiTia) (12) 

z = (L^Tl 2 ) (13) 

and 2^ — 3 invariants: 

X a = L1.L2 (14) 

Xij = Li.Lj (15) 

X 2j = L 2 .L 3 (16) 



for j ^ 3. We also showed how the power spectrum and 
bispectrum can be written in terms of the Xij and A. 

Note that the Multipole Vectors are "headless" vectors: 
they do not have a defined direction. This can be clearly 
seen in Eg uation 1101 where any change in the sign of a vec- 
tor can be absorbed into the scalar A. Therefore, we have 
to be careful in the definition of our orthonormal frame. We 
impose possible directions on the "anchor vectors" by insist- 
ing that their directions satisfy X12 = L\.L 2 > 0. Therefore, 
these vectors can still point in either direction but their di- 
rections are not independent, and there are only 2 possible 
combinations. The z-axis is now uniquely defined, and the 
x, and y-axes are up to a rotation by tt around the z-axis. 
We simply reduce the allowed range for the Euler angles of 
the axes. 

We now have an orthonormal frame defined by our an- 
chor vectors; we call this the Multipole Frame. The remain- 
ing vectors' information is contained in their angles. The 
minimum number of dot products needed to contain all the 
information about the vector positions is the 2^ — 3 invari- 
ants Xij . So we see that information about m-preference and 
Gaussianity are contained in the angles between vectors, i.e. 
their relative spread. 

Gaussianity and SI determine the probability distribu- 
tion of these angles (see lDennisl d200 4b) for an analytical ex- 
pression). Departures from the expected spread would indi- 
cate a departure from SI and/or Gaussianity. It is not possi- 
ble to probe the issues of m-preference and non- Gaussianity 
separately, however in some extreme cases of m-preference 
the vectors form clear structures. For pure m-preference 
there is a frame where [i — m) of the vectors will line along 
the z-axis of this frame, and the remaining m vectors will 
form a 2m regular polygon in the x-y plane. We refer to this 
as a "Handle and Disc" structure, with the aligning vectors 
making the "Handle", and the planar vectors forming the 
"Disc". 

Note that this method involves all the information con- 
tained in the multipoles, as opposed to methods that explore 
just one statistic e.g., But- This obviously has many advan- 
tages as it allows for a much more thorough investigation. 
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Figure 1. The Euler angles found for the Multipole frames of 1000 simulations, for multipoles 1 = 2- 
sky, and the dashed line is with the KpO masked applied. 



20. The solid line is for the whole 



3.3 Anchors 

We now explain further how to select the anchor vectors 
essential for our separation of SI and non-Gaussian issues. 
We need an internal way of selecting these, i.e., a method 
that does not depend on any outside coordinate frame. How- 
ever, we want the process to be repeatable so selecting them 
at random (as we've done before) is not satisfactory. For 
instance we may want to compare results for different fre- 
quency channels, so that the algorithm must be unambigu- 
ous as well as invertible. 

One such way would be to define the vectors as the 
two most aligned, or most orthogonal. This would select a 
pair, but would not provide them or the remaining vectors 
with an order. This is not a problem in defining the frame. 
However, with no ordering of the Xij, X2j we cannot assess 
the agreement of individual values with the simulations; we 
can only study them as a set. Ideally we want to use every 
piece of information, and therefore we look for a way to order 
the vectors. 

We achieve this in the following way. From the set of £ 
unordered vectors {vi}, for every vector we sum the square 
of its dot products with all other vectors: 

Ki = J2(X^f (17) 

3 

where Xij = Vi.Vj. Every vector now has a corresponding 
value of Ki, and we can order them. For example let L\ be 
the vector with the maximum Ki value, and L2 equal the 
vector with the second highest Ki value, etc... so Li is the 
vector with the lowest Ki value. We can now calculate the 
2^ — 3 Xij values, and use the vectors L\ and L2 to define 
the frame. 



4 A SIMPLE APPLICATION 

A sky cut {e.g., imposed by the KpO mask) induces 
Anisotropy onto the sky. To investigate the effectiveness of 
our method we look at how well this Anisotropy is detected. 
We simulate 1000 Gaussian skies with noise and beam, and 
we look at the orthonormal frames before and after we ap- 
ply the KpO mask. We find the multipole frame {x, y,z}, 
as defined above, for multipoles £ — 2 — 20, and record the 



corresponding Euler angles {<f>,9,if)) 1 . As mentioned in Sec- 
tion 13.21 due to the Multipole Vectors having no defined 
direction, the definition of the orthonormal axes is defined 
up to (x,y) — > (— x, — y). Therefore there is an equiva- 
lence between orthonormal axes {x, y,z} ~ {— x, — y,z}, 
which translates into an equivalence between Euler angles 
(4>,9,ip) ~ (4> + n,9,ip). Thus the allowed ranges for the 
Euler angles of our Multipole Frames are cf> g [0, 7r],0 £ 
[0,7r],V G [-7T,7r]. 

In Figure Q we plot histograms of the distributions 
found for the Euler angles. We see that without any sky- 
cut the distributions are uniform as expected. However, we 
also clearly see how the sky-cut then does induce anisotropy. 
The 4> and cos(9) distributions become non-uniform: 9 peaks 
around 7r/2 (2-axis in galactic plane), and (j> also around 
7r/2 (rr-axis aligns with galactic poles). This is unlike other 
measures, such as the bispectrum, which doesn't get largely 
effected by sky cuts. 

In fact, one can predict this result. In terms of the 
Fourier expansion in spherical harmonics, the effect of the 
mask is to remove power from the m ~ ±£ in galactic coordi- 
nates. Therefore we are equivalently inducing m-preference 
towards the lower m's. As discussed in Section 13.21 pure 
m-preference leads to a particular "Handle and Disc" con- 
figuration for the vectors, with £ — m vectors in the handle 
and m in the disc. Therefore, for this induced preference 
towards the lower m's we would expect a handle and disc 
structure (if only a subtle one) with most of the vectors in 
the handle. Our method of selecting the anchor vectors will 
then favour those in this handle as they will maximise Ki 
through their small inter-angles with other handle members. 
If the anchor vectors belong to the handle, then the a;-axis 
will align with this handle and therefore point to the galac- 
tic poles, the t/-axis will not be well defined, and the 2-axis 
will be perpendicular to the handle and therefore lying in 
the galactic plane. 



1 We are using the "x-convention" , where the Euler angles 
(</>, 0, ip) represent first a rotation by <j> around the z-axis, 
then by 8 around the x-axis, then by if) around the 2-axis 
(again). Further explanation can be found at MathWorld: 
http: / /mathworld. wolfram, com 
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1 = 2 1 = 3 £ = 4 




£ = 5 i=6 1 = 7 




£ = 8 £ = 9 £ = 10 

Figure 2. The Multipole Vectors for multipoles £ = 2 — 10 of the KpO masked Inverse-Noise-Squared coadded WMAP map, coadd. 
The vectors are headless, so each vector is represented by two white dots, where it intersects with the unit sphere. 
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£=5 £=6 £=7 
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5 RESULTS 

5.1 The datasets and methodology 

We investigate the multipole vectors of the CMB using the 
WMAP data. We look at different representations of the 
data in the form of four different maps: 

1. An Inverse-Noise-Squared Coadded map (coadd). We 
take the 8 foreground cleaned maps in the Q, V, and W 
bands, provided by WMAP 2 . We combine these 8 data ar- 
rays (DAi) using coefficients that maximally reduce the 
noise in the final map: 

T(n) = ^„ /; (18) 

where Ti(n) is the sky map for the DA i with the foreground 
galactic signal subtracted, and 

Npix 

oi = v^ ? w (19) 

71 = 1 

af(n) = ao,i/Nobsi(n) (20) 

(Tpi is t ne noise per observatio n for DA i, whose values 
are given by iBennett fc all iFoiKii. We apply the KpO mask 
to this map feennett fc alll2003l) to cut possible remaining 
galactic contamination near the disk, and then we remove 
any residual monopole. 

2. The Inter-Linea r-Combination map (ILC), produced by 
the WMAP team feennett fc aJl2003) . This map is formed 
from a combination of smoothed maps, with weights cho- 
sen to maintain CMB anisotropy signal while minimizing 
the Galactic foreground. The WMAP team advise that the 
ILC map be used as a visual tool, and not used for CMB 
studies. We include it here for interest, as other groups 
have claimed interesting detection s of anomalies in this 
map(|Cogi^^uterer^^aAma^|200^ ; _ 

3+4. Tfegmeulc^Odlveha^C^ste^Ianiiltonl (I2003T) produced 
their own foreground cleaned maps 3 . Their approach did 
not assume a specific power spectrum, and these maps con- 
tain less foreground signals and noise than the ILC map. 
For the lowest multipol es the maps are clean enough that no 
galac tic cut is required llTegmark. Oliveira-Costa. Hamiltonl 
2003). Our third map is their cleaned map (TOHcl), and 
our fourth map is their Wiener filtered map (TOHw). 

We will compare findings to those from simulations. Our 
simulations are of Gaussian fluctuations, with the noise and 
beam properties of the coadd map. We perform 1000 of 
these simulations with the KpO mask, and 1000 without. In 
the case of the coadd map we always compare the results 
to those with the KpO mask added, and for the other ILC, 
TOHcl, and TOHw maps we compare their results to those 
without the mask. We have only done simulations with the 
noise and beam properties of the coadd map to simplify 
matters. We find this acceptable as at the low-£ multipoles 
we are concentrating on, the difference in the levels of noise 
and resolution seen in the four maps are insignificant. 

Throughout our analysis we use software produced 
by the HEALPix collaboration , iGorski. Hivon. Wa ndclt 

2 Available at |http: / /lambda. gsfc. nasa.gov/] 

3 Available at www.hep.upenn.edu/ max/wmap.html 



(1999) 4 . To compute the Multipole Vectors we use codes 
produced bv lCopi. Huterer. Starkmanl i2003tl 5 . 



5.2 The Multipole vectors 

In Figures I2l.il we display in Mollweide projection the 
multipole vectors for map coadd and ILC, for multi- 
poles £ = 2 — 10 (for these multipoles the vectors for 
maps TOHcl and TOHw by eye look very similar to 
the ILC map results). Recall that these vectors are head- 
less, so we show dots at both ends. In Figure [3] one can 
see the planarity of the £ = 3 multipole vectors, and 
the alignment of this plane with the galactic plane and 
the plane of the £ = 2 mul ti pole vectors, as reported 
bv | 1i v eira-Costa fc al| (|2003}l: [Cpp^ Huterer. Starkmanl 
j2003ft : ISchwarz fc alH2003ft : IWeeksl (l2004h . However, a word 
of warning, this visual projection can be misleading as other 
significant features would not be so apparent by eye. In this 
projection our eye easily picks out any planarity in the galac- 
tic plane. 

We now use our prescription of defining "anchor vec- 
tors" to investigate the multipole-frames and the multipole- 
dots for each map. Therefore the analysis of CMB data now 
splits into two areas: 

1. The Multipole Frames. We search for inter-^ correlations 
by looking for correlations between these frames. SI demands 
that the frames be uniformly orientated (apart from the ef- 
fect of sky cuts, noise, etc). Any correlation between frames 
of different multipoles, apart from that seen in the simu- 
lations, would indicate statistical anisotropy. The simplest 
case would be an alignment of the frames - clearly indicating 
one preferred frame. We investigate this issue for our four 
maps in Section [5.31 

2. The Multipole Invariants. For each multipole we extract 
the 2£ ~ 3 invariants, Xij. From simulations we can extract 
their expected distributions assuming SI and Gaussianity 
taking into account additional factors such as noise and 
beam. We investigate this in Section f5. 41 



5.3 The Multipole Frames 

For each of the four maps we do the following. For each mul- 
tipole in the range £ — 2 — 20 we find the Multipole Vectors, 
and order them as described in Section 13.31 We define L\ 
and L2 as the anchor vectors, and use these anchor vectors 
to define an orthonormal frame as described in Section 13. 21 
Under the hypothesis of Statistical Isotropy we expect 
these multipole frames to be uniformly distributed (except 
for the anisotropy induced by the galactic cut in the case 
of the coadd map). One cannot assess how uniformly dis- 
tributed a frame is one at a time as clearly every orientation 
is just as likely and so no significance can be assigned to any 
particular result. One can only ask how probable it is that 
a number of frames came from a uniform distribution. To 
do this we use a Kolmogorov-Smirnov test. This provides a 
goodness of fit test for a statistical distribution. It looks at 



4 Available at http://www.eso.org/science/healpix 

5 Available at www.phys.cwru.edu/projects/mpvectors/ 
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cos(6) 




l = 2 - 10 


coadd 
ILC 
TOHcl 
TOHw 


0.335(9.0) 
0.079(92.3) 
0.370(5.1) 
0.259(25.2) 


0.205(51.5) 
0.256(35.5) 
0.174(73.8) 
0.145(87.7) 


0.360(12.0) 
0.284(26.2) 
0.138(90.3) 
0.249(40.5) 


I = 2 - 20 


coadd 
ILC 
TOHcl 
TOHw 


0.213(7.8) 
0.130(49.0) 
0.133(44.5) 
0.133(44.5) 


0.187(25.3) 
0.070(98.4) 
0.225(14.8) 
0.075(94.1) 


0.219(17.4) 
0.069(95.9) 
0.109(80.6) 
0.136(57.1) 



Table 1. D-statistics from cumulative frequency distributions of 
the Multipole Frame Euler angles from multipoles I = 2 — 10 and 
I = 2 - 20. 



the discrepancy between the cumulative frequency of obser- 
vations and the expected cumulative frequency. The maxi- 
mum discrepancy is called the D-statistic. For large enough 
sample sizes there is an analytical expression for the critical 
values of the D-statistic, over which a null hypothesis should 
be rejected. 

We calculate the cumulative frequency over a range of 
multipoles of the Euler angle observations, <f>,cos(9),ip. We 
do this for two multipole ranges, £ = 2 — 10 and £ = 2 — 20. 
We compare the distributions we find for the maps to those 
from simulations. As expected, the distributions the simula- 
tions return are near uniform. For example, in Figure 0] we 
visualise the cumulative frequencies of the coadd map re- 
sults, accumulated over multipoles £ = 2 — 20. The expected 
distribution returned by the simulations is also shown. The 
D-statistic is the maximum discrepancy between the two 
lines. 

For each of the four maps, and for both the £ ranges 
we find the maximum discrepancy, the D-statistic, that the 
3 Euler angles find. As our sample sizes are relatively small 
(9 observations for £ — 2 — 10 range, and 19 for £ = 2 — 20) 
we do not use analytical expressions for the critical values of 
the D-statistic: we assess the significance of the D-statistics 
we observe by comparing them to all those the simulations 
found. In TableQwe list the D-statistics found and, in brack- 
ets, the percentage of simulations that found larger values. A 
significant departure from the null hypothesis would result 



in a high D-statistic, and therefore a low percentage value. 
We find no significant evidence for any departure from the 
expected uniform (nearly uniform in the case of the coadd 
map) distribution for the multipole frames. Therefore, using 
our Multipole Frames we find no evidence for a departure 
from SI for multipoles I = 2 - 20. 

5.4 The Multipole Invariants 

We now turn our attention to the invariant degrees of free- 
dom. As discussed in Section r3.ll the invariants of one multi- 
pole will probe Gaussianity, and the m-preference associated 
with SI. In terms of the Multipole Vectors this information 
is stored in t heir rela t ive spa cing - the angles between one 
and another. IPennisI (l2004bl) computed the analytical ex- 
pression for the expectation of the general dot product for 
a given multipole. 

For each multipole we will not look at all the dot prod- 
ucts between all the vectors as this contains redundant in- 
formation. We will use our prescription of defining anchor 
vectors, and observe the 2^ — 3 invariants Xy , as described 
in Section |3~T1 

For each map, and for each multipole £ = 2 — 20 we 
calculate the 2£ — 3 invariants Xij . We compare their values 
to those from simulations. In Figure|S]we display the results 
found for the invariants X12, X13, and X23. We also show 
the Xij ± a values that the simulations returned. 

In these plots we can observe potentially anomalous fea- 
tures. For example in the Xn plot the £ = 5 value is very 
low for the ILC, TOHcl, and TOHw maps, and the £ = 13 
value is low for the coadd map. The A13 plot again shows 
out-lier values for £ — 5 in the ILC map, and £ = 13 for the 
coadd map. 

To assess the significance of the outliers, we compare 
the values of all the invariants Xij for multipoles £ — 2 — 
20 to those from simulations: for each value we find the 
number of simulations that find a lower value. We highlight 
as anomalous any result where either >99% or <1% of the 
simulations found lower values, and we record just these 
results in Table [5] 

Before we take any of these results too seriously, we 
must note that a certain of number of anomalous values are 
expected. That is the same with any statistic. For each multi- 
pole there are 2^ — 3 invariants that we measure. Therefore, 
when we look at the multipoles £ — 2 — 20 for each map 
we observe a total of 361 invariants. We have highlighted 
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Figure 5. The X12, X13, and X23 results for the coadd map (bold line), the ILC map (short-dash), the TOHcl map (long-dash), and 
the TOHw map (dash-dot). Also shown is the mean, ±1<t for the simulations (dotted), with and without the KpO mask added. 
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Table 2. Anomalous multipole invariant values, at the 98% signif- 
icance level. The percentage of simulations that find lower values. 



the values outside the 98% confidence interval, and so we 
would expect 2% of our observations to be highlighted: ~7 
results for every map. We see in the tables that the coadd 
map finds 8 "anomalies", the ILC finds 9, TOHcl finds 6, 
and TOHw finds 9. We also expect to find these ~7 values 
in the higher multipoles, as these return more invariants. 

Bearing this in mind, we highlight as particularly inter- 
esting the "anomalies" found in the low multipole £ = 3, for 
maps TOHcl and TOHw. This result shows a particularly 
low dot product between the second and third vectors, and 
therefore an unusually large angle. 

Also of significance are the multipoles where 3 "anoma- 
lous" values are found: £ = 13 for the coadd map, and 
I = 18 for the ILC map. 

However, we conclude that through analysis of the 
Multipole Vectors there is no evidence for inherent non- 
Gaussianity or m-preference in the multipoles 1 = 2 — 20. 



6 DISCUSSION 

With improving CMB observations, the job of analysing the 
data also requires improved methods. It is unsatisfactory to 
investigate numerous statistics that bear no intuitive inves- 
tigative properties. In particular, the issues of SI and Gaus- 
sianity should be probed separately. 

In this paper we ha ve applied a method proposed 
in lLand. Magueiiol (|2004rJl to the WMAP data. The method 
uses all the information in a map, and separates out the 
issues of SI and Gaussianity by defining an orthonormal 
frame (Multipole Frame) for each multipole, and a set of 
invariants (Multipole Invariants). In theory, the robustness 
of this method makes it ideal for a thorough investigation 
of anomalies in the CMB. 

We analysed the Multipole Vectors of four different 
CMB maps, all derived from the WMAP data. We looked at 
the multipoles I — 2 — 20, and for each multipole we found 
a Multipole Frame and the Multipole Invariants. To inves- 
tigate the Statistical Isotropy of the data we analysed the 
orientation of the Multipole Frames over 2 separate multi- 
pole ranges, I = 2 — 10, and £ = 2 — 20. We found no evidence 
that the Euler angles of these frames significantly favoured 
any particular values. Therefore we found no evidence for 
a departure from SI. To investigate the Non- Gaussianity of 
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the data we analysed all the 361 Multipole Invariants for 
multipoles t — 2 — 20 and compared their values to those 
from simulations. We recorded the extreme values, and saw 
that each map found approximately the expected amount. 
Therefore we found no evidence for Non-Gaussi anity. 

In light of previously detected anomalies iCole s&^ajl 
20031: ICopi. Huterer Starkmanlliool lLand. Magueiidl200a 
Oliveira-Costa fc all2003t iRalston. JairJl2004 ISchwarz fcal 
20031 : iTegmark. Oliveira-Costa. HamiltorJl2003r) . the failure 
of our method to return any significant results raises a cou- 
ple of interesting points. The effectiveness of our method 
comes into question, and in practice we have seen that our 
method is victim to discontinuous noise. The method pivots 
around the selection of Anchor Vectors. Small fluctuations 
in the positions of the Multipole Vectors alters their order- 
ing and therefore has a large effect on our method. The 
Anchor vectors can then differ, and this induces a discontin- 
uous noise in the Euler angles of the Multipole Frames and 
in the Multipole Invariants. Also, different ways of order- 
ing the vectors will make the method effective at detecting 
different features. 

Another significant issue raised is that of priors. In this 
paper we limit ourselves to the multipoles I — 2 — 20, but 
we do not focus any further or assume priors about any 
range. Therefore we may overlook interesting features in a 
few multipoles because their anomalous behaviour will be 
diluted by the well behaved values from other multipoles. In 
only focusing on interesting multipoles one has assumed a 
prior. In this paper we have applied our method to look for 
more overall features. 

We conclude that the Multipole Vector method with 
the anchor vectors is technically ideal, but in practice is 
very limited by noise. In its application we have found no 
evidence for anisotropy or non-gaussianity. However, we feel 
that in light of other reported results this is because our 
method overlooks subtle features in the data. What we gain 
in thoroughness, we loose in sensitivity. We find no evidence 
for inherent anisotropy or Non-Gaussianity. 
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